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A variety of structural abnormalities have been described in post traumatic stress disorder ( PTSD), but only a few 
studies have focused on cortical thickness alterations in recent onset PTSD. In this study, we adopted surface- 
based morphometry (SBM), which enables an exploration of global structural changes throughout the brain, in 
order to compare cortical thickness alterations in recent onset PTSD patients, trauma-exposed subjects but 
without PTSD, and normal controls. Moreover, we used region of interest (ROI) partial correlation analysis to 
evaluate the correlation among PTSD symptom severity and significant changes of cortical thickness. The wide- 
spread cortical thickness reduction relative to the normal controls were found in bilateral inferior and superior 
parietal lobes, frontal lobes, hippocampus, cingulate cortex, and right lateral occipital lobes in trauma survivors, 
whereas cortical thickness was only increased in left calcarine cortex in PTSD group. The average cortical thickness 
of hippocampus and cingulate cortex decreased by 10.75% and 9.09% in PTSD, 3.48% and 2.86% in non PTSD. We 
further demonstrated that the cortical thicknesses of bilateral ACC and PCC, superior frontal lobes, and hippocam- 
pus are negatively correlated with CAPS scores in all trauma survivors. Our study results suggest that stress 
widens cortical thinning regions and causes more serious effect in recent onset PTSD than non PTSD. It also 
shows that the cortical thinning in recent onset PTSD predicts the symptom severity. 

© 2013 The Authors. Published by Elsevier Inc. All rights reserved. 



1. Introduction 

Neuroimaging studies have identified a number of functional and 
structural alterations in recent onset post-traumatic stress disorder 
(PTSD). Structural neuroimaging studies about recent onset PTSD 
focused primarily on gray matter volume alterations based on voxel 
based morphometry (VBM) and manual based hippocampal and 
amygdale volumetry. These studies have shown that recent onset 
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PTSD is associated with smaller gray matter volume of limbic structures 
(Chen et al., 2006; Corbo et al., 2005; Wignall et al., 2004; Zhang et al., 
2011). However, there are inconsistencies in previous study results 
(Bonne et al., 2001; Jatzko et al., 2006). Compared to the normal 
controls, decreased gray matter volumes or densities were identified 
in left or right hippocampus, left parahippocampal gyms, bilateral 
calcarine cortex, left or right anterior cingulate cortex (ACC), bilateral 
or left insular cortex in recent onset PTSD subjects (Chen et al., 2006; 
Corbo et al., 2005; Zhang et al., 2011). Volumetric studies identified 
the smaller right hippocampus volume and whole brain volume 
(Wignall et al., 2004), while the longitudinal volumetric study did not 
detect any hippocampus volume changes (Bonne et al., 2003). As men- 
tioned in our previous study, numerous confounding factors can be 
possible reasons for the difference of VBM study results, such as a differ- 
ent study method or trauma type or trauma duration time, or subtle 
cerebral cortex impairment (Landre et al., 2010; Zhang et al., 201 1 ). 

Cortical thickness is a reflection of the size, density, and arrangement 
of neurons, glial cells and nerve fibers (Narr et al., 2005). It is a parameter 
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relatively invariant to brain size during mammalian evolution, as well as 
the content of, the neuropile (Keep and Jones, 1990). Although cellular 
characteristics cannot be quantified directly in imaging data, cortical 
thickness was thought to reflect cytoarchitectural abnormalities more 
closely than cortical volume does (Keep and Jones, 1990; Thompson 
etal., 2003). 

VBM analyses are particularly sensitive to the degree of smoothing, 
differences in registration, and choice of normalization template 
(Bookstein, 2001 ; Jones et al., 2005; Park et al., 2004). Thus, the sur- 
face-based morphometry (SBM) analysis approach has been put for- 
ward as an alternative method for probing into the cortical gray 
matter changes, which also allows for the contributions of gray 
matter thickness and regional surface area to be defined independently. 
The FreeSurfer software is a feasible method to precisely quantify and 
characterize the cortical thickness in different brain regions (Fischl and 
Dale, 2000; Makris et al., 2006; van der Kouwe et al., 2008), and has 
been used in PTSD studies (Dickie et al, 2012; Geuze et al., 2008; 
Hunter et al, 2011; Kuhn et al, 2011; Landre et al, 2010; Liu et al, 
2012; Lyoo et al, 201 1 ; Woodward et al, 2009a). 

In previous studies, imaging data support the neurocircuitry model 
of PTSD that emphasizes the functional relational relationship between 
the amygdale, ACC and hippocampus. However, most of previous stud- 
ies focused on chronic PTSD and non PTSD cortical changes after trauma, 
and cortical thinning were found in parahippocampal gyrus, superior 
temporal cortex, and lateral orbital frontal cortex (Woodward et al, 
2009b). ACC was a steady predication for PTSD recovery (Dickie et al, 
2012), and cortical thickening after trauma exposure (Lyoo et al, 
2011). However, there are many confounding factors may influence 
the results of chronic PTSD, such as substance use/abuse, trauma 
exposure time and trauma type. To reveal the subtle alteration in recent 
onset PTSD, trauma survivors from a single trauma event, with homoge- 
neity in demographic variables and traumatic type, intensity, and dura- 
tion of exposure, and avoided most comorbidity factors, may offer a 
distinct advantage in examining psychological consequences of PTSD 
and stress on brain structure. 

The purpose of this study is to detect the pathogenesis of recent 
onset PTSD and stress by utilizing the SBM method. Based on the results 
of previous studies, we hypothesize that there may be explicit regional 
cortical thickness differences between the recent onset PTSD and 
trauma exposed non PTSD as compared with the healthy controls, and 
that the regional cortical thickness abnormalities may be associated 
with the severity of PTSD symptom. 

2. Methods 

2.3. Subjects 

On July 29, 2007, a severe coalmine-flood disaster occurred at 
Zhijian Coalmine in Shanxian County in Henan Province. Sixty-nine 
male miners were rescued after a 75 hour ordeal in the darkness. 
Three months after the disaster, we attempted to contact all the 69 
miners rescued. Forty-eight of them (69.6%) were successfully found 
and were willing to participate in the study. Six months later, 17 of 
these 48 coal mine disaster survivors met DSM-IV (Anon, 1994) diag- 
nostic criteria for PTSD. 10 survivors and 5 survivors' wives with PTSD, 
10 survivors without PTSD, and 25 normal controls were enrolled in 
the MR1 study. The study protocol was approved by the Medical Ethical 
Committee of Xijing Hospital of the Fourth Military Medical University. 
Written voluntary and informed consent was obtained from all partici- 
pants before the study. Each psychiatrist was trained to expertly use 
the DSM-IV and the Structured Clinical Interview for DSM-IV (SCID) 
(Spitzer et al, 1996) to diagnose PTSD. None of the subjects met diag- 
nostic criteria for major depression, schizophrenia or bipolar disorder. 
The identified PTSD subjects had never received psychiatric treatment 
for their condition. Moreover, none of the subjects had a history of treat- 
ment with psychotropic drugs or of substance abuse. All subjects were 



free of metallic implants, neurological disorders and major medical con- 
ditions. The severity of their symptoms was assessed with the Clinician- 
Administered PTSD Scale (CAPS). The elapsed time between the trau- 
matic event and the MRI scan ranged from 187 to 190 days. 

2.2. Magnetic resonance imaging acquisition 

All MRI scans were performed on a 3.0 T MR scanner (MAGNETOM 
Trio, Siemens AG, Erlangen, Germany) located in the Department of 
Radiology, Xijing Hospital. A high-resolution three-dimensional magne- 
tization prepared rapid acquisition gradient echo (MPRAGE) Tl- 
weighted sequence was used to acquire MR images of the whole 
brain (176 sagittal slices). The other MR imaging parameters applied 
in this study were as follows: TR = 1900 ms, TE = 2.26 ms, TI = 
900 ms, flip angle = 9°, acquisition matrix = 256 x 256, field of 
view = 220 x 220 mm 2 and 1.00 mm slice thickness with no inter- 
slice gap and isotropic resolution 1.0 x 1.0 x 1.0 mm 3 . 

2.3. MR image processing 

Cortical reconstruction and parcellation were performed with the 
FreeSurfer V4 package (v4.0.2) image analysis suite (Fischl and Dale, 
2000), which is documented and freely available for download online 
(http://surfer.nmr.mgh.hai-vard.edu/). This method is fully described 
(Du et al, 2007) and has been validated via both histological and man- 
ual measurements (Bischoff-Grethe et al, 2007; Pantazis et al, 2010). 
For cortical analysis, the gray matter-white matter interface of each 
hemisphere was segmented, tessellated, corrected for topological errors 
and inflated in order to unfold the cortical surface. Each individual 
hemispheric surface was then registered to a surface-based template 
using its sulco-gyral pattern to allow inter-subject analysis (Fischl and 
Dale, 2000). Cortical thickness was calculated from both intensity gradi- 
ents across tissue class and continuity information, as the closest dis- 
tance between the gray/white matter boundary and the gray/cerebro 
spinal fluid interface at each vertex of the tessellated surface (Fischl 
and Dale, 2000). Automated topology correction was performed as nec- 
essary (Fischl et al, 2001 ). We inspected each dataset at multiple points 
within the processing stream. We verified the Talairach transform, the 
accuracy of the skull strip, the accuracy of the white matter and pial 
surface segmentation, and finally inspected the cortical surface recon- 
struction for topological defects. Error correction was performed as 
necessary. Cortical thickness was measured as the distance between 
the gray/white matter boundary and the pial surface at each point on 
the cortical mantle. The data for each participant was resampled to an 
average participant and surface smoothing was performed using the 
'qcache' command and a 10 mm full-width half-maximum Gaussian 
kernel, prior to statistical analysis. In this study, we get the cortical 
thickness values from the statistical output of FreeSurfer cortical 
parcellation, called "aparc.a2005s. stats", which uses the Desikan- 
Killiany atlas, and has 75 regions. 

Once the processing stream was satisfactorily completed, the QDEC 
(Query, Design, Estimate, Contrast) which is a FreeSurfer application 
intended to aid researchers in performing group averaging and inference 
on the morphometry data (cortical surface and volume), was used to 
explore cortical thickness differences between groups and the relation 
of cortical thickness to the cognitive performance variables(CAPS scores). 
Statistical parametric maps of the entire cortical mantle were generated 
to show group differences. 

2.4. Statistical analysis 

2.4.1. Cortical thickness analysis 

A general linear model (GLM) was estimated at each vertex across 
the cortical surface, with cortical thickness as dependent variable, diag- 
nosis and sex as categorical predictors, and age as continuous predictor. 
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This procedure allows for a generation of statistical parametric maps to 
be plotted. 

The maps show the distribution of p values for pair-wise compari- 
sons between the diagnostic categories, as defined by the following 
contrasts: (1) PTSD versus the normal controls; (2) PTSD versus non- 
PTSD subjects; (3) non-PTSD subjects versus the normal controls. To 
correct multiple comparisons, p-maps were subject to an expected 
False Discovery Rate (FDR) correction of p < 0.05. 

2.5. Correlation analysis 

2.5.2. Correlations between ROI thickness and CAPS 

The mean cortical thickness value of regions with significant cortical 
thinning was extracted from all trauma survivors. The mean regional 
cortical thickness and CAPS values were imported into the SPSS17.0 
software (SPSS Inc., Chicago, IL) for correlation analysis, with age and 
gender being taken as covariate. 

2.5.2. Cortical thickness change rates 

To reduce the impact of individual differences, we measured the 
PTSD cortical thickness change rate in both PTSD and non PTSD, as 
cortical thickness change rate = (mean normal control — mean 
PTSD) x 100 / normal control cortical thickness. The same equa- 
tion was also used for non PTSD cortical thickness change rate. 

3. Results 

3.1. Demographics 

All subjects are right-handed, and there are no significant differences 
between groups in gender and age. All subjects in this study came from 
the same community and did not differ significantly in socioeconomic 
status (Table 1). As expected, the PTSD group had higher scores on the 
CAPS than the non-PTSD group (PTSD group: 84.25 ± 14.29; non- 
PTSD group: 33.15 ± 20.12) [F = 34.75; df = 23, p = 0.000]. 

3.2. Cortical thickness 

Global thickness measurements were examined by pair-wise com- 
parisons, with age and gender as independent variables. There were 
no significant differences between groups in terms of whole brain aver- 
age cortical thickness. But the cortex was significantly (FDR< 0.05) 
altered both in PTSD and non-PTSD, compared to the healthy controls 
(Table 2). 

3.3. Cortical thickness analysis in groups 

3.3.3. PTSD vs. healthy control 

3.3.1.1. PTSD < healthy control. Cortical thickness was significantly 
decreased in PTSD compared with the healthy control subjects in 

Table 1 

Demographic characteristics of recent onset PTSD, non PTSD and normal controls. 



Variable PTSD Non-PTSD Control 

(mean + SD) (mean + SD) (mean + SD) 



Age 


41.6 ± 5.93 


34.3 ± 7.37 


36.6 ± 3.25 


Sex (male/female) 


10/5 


10/0 


22/3 


Marital status (Ma/Unma) 


14/1 


10/0 


23/2 


Smoked in previous 1 year 


11/4 


8/2 


10/15 


The length of being a miner 


4.7 ± 1.3 


6.2 ± 1.1 




The length of education 


6.67 ± 2.06 


7.25 ± 2.38 


17 ± 2.14 


BMI 


26.34 ± 1.93 


25.47 ± 1.55 


25.12 ± 1.23 


CAPS 


84.25 ± 14.29 


33.15 ± 20.12 




Global cortical thickness 


2.713 ± 0.109 


2.721 ± 0.133 


2.726 ± 0.193 



PTSD, post-traumatic stress disorder; SD, standard deviation; CAPS, Clinician-Administered 
PTSD Scale; BMI, Body Mass Index; Ma, married; Unma, unmarried. 



multiple regions across the cerebral cortex. In the left hemisphere, 
the cortical thinning regions were in the superior frontal gyrus, hip- 
pocampal gyms, inferior parietal lobes, anterior cingulate cortex and 
posterior cingulate cortex. In the right hemisphere, the thinning regions 
were in the superior frontal gyrus, rostral middle frontal gyrus, hippo- 
campal gyrus, superior parietal lobes, inferior parietal lobes, lateral 
occipital gyrus, anterior cingulate cortex, and posterior cingulate cortex 
(Fig. 1). 

3.3.1.2. PTSD > healthy control 

There was cortical thickening only in the left calcarine cortex ( Fig. 1 ) . 

3.3.2. PTSD vs. non-PTSD 

3.3.2.1. PTSD < non-PTSD. In the left hemisphere, the cortical thinning 
regions were in the superior frontal lobes, hippocampus, the ACC and 
PCC. In the right hemisphere, the cortical thinning regions were in the 
lateral occipital gyrus, rostral middle frontal gyrus, hippocampus, the 
ACC and PCC. 

3.3.2.2. PTSD > non-PTSD 

There was cortical thickening only in the left calcarine cortex (Fig. 2). 

3.3.3. Non-PTSD vs. healthy control 

The cortical thickness did not show significant difference between 
the two groups. There were no significant findings for this contrast. 

3.3.4. ROI analysis 

3.3.4.1. Correlations between regional cortical thickness and CAPS. There 
were negative correlations of cortical thickness with CAPS in the 
cingulated cortex, inferior parietal cortex, superior frontal cortex 
and hippocampus (as shown in Table 3, Fig. 3). 

In ROI analysis, compared with the healthy controls, the average 
cortical thickness of the hippocampus and the cingulate cortex de- 
creased by 10.75% and 9.09% in PTSD, by 3.48% and 2.86% in non 
PTSD. The detail showed that the left part of the hippocampus de- 
creased by 10.13% in PTSD and by 3.61% in non PTSD; the right part 
of the hippocampus decreased by 11.34% in PTSD and by 3.35% in 
non PTSD; the left part of the anterior cingulate cortex decreased 
by 11.76% in PTSD and by 3.89% in non PTSD; the right part of the ante- 
rior cingulate cortex decreased by 11.57% in PTSD and by 3.03% in non 
PTSD; the left part of posterior cingulate cortex decreased by 6.11% in 
PTSD and 1.51% in non PTSD; the right part of the posterior cingulate 
cortex decreased by 6.95% in PTSD and by 3.03% in non PTSD (and the 
details list in Table 2). 

4. Discussion 

The main findings of the present study are: 1 ) There is widespread 
regional cortical thinning in recent onset PTSD patients relative to non 
PTSD subjects and the normal controls. The region of cortical thinning 
in recent PTSD patients is much larger than that of non PTSD subjects. 
The cortical thinning region in PTSD is primarily in the bilateral hippo- 
campus, ACC and PCC, left inferior and right superior parietal lobes, 
left superior frontal lobe, and right superior rostral middle frontal 
gyms. The cortical thickening in recent onset PTSD patients is only in 
the left calcarine cortex. 2) There is no significant cortical thickness 
difference between non-PTSD and the normal control. The mean global 
cortical thickness of the three groups did not show significant differ- 
ence. 3) In ROI analysis, the cortical thickness of the bilateral hippocam- 
pus, the cingulate cortex and frontal gyrus negatively correlated with 
the CAPS score in all trauma victims. 4) The cortical thinning rate of 
the bilateral hippocampus and cingulate cortex in recent onset PTSD pa- 
tients is much larger than that in non-PTSD subjects. These results imply 
that severe stress has profoundly selective effects on the medial-frontal 
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Table 2 

Cortical thickness in recent onset PTSD, non-PTSD and normal control. 



Cortical area 


Side 


Mean (SD) in mm 






PTSD vs. non-PTSD 


PTSD vs. Healthy 


PTSD 

In — 1^1 

(_n — lb) 


1NOI1-1 1 jU 
[Vl — IUJ 


neaiLiiy 

tn — 9^ 

l n — zdj 


T value 


P value 


T value 


P value 


Cortical thickness decreasing (PTSD < 


control) 
















Superior frontal gyrus 


Left 


2.389 (0.128) 


2.422 (0.125) 


2.609 (0.251) 


— J.OlJ 


u.uzz 




n fine 




Right 


2.411 (0.247) 


2.477 (0.148) 


2.543 (0.268) 


— z.ooy 


U.UJO 


A 1(W 


U.Ulj 


Rostral middle frontal gyrus 


Left 


2.147 (0.220) 


2.157 (0.196) 


2.252 (0.181) 


1 P.AA 




1 oa 


IN j 




Right 


2.199 (0.217) 


2.228 (0.234) 


2.331 (0.195) 


— z.zyjj 


U.U4Z 


— 


U.UjZ 


Hippocampal gyrus 


Left 


2.291 (0.208) 


2.457 (0.215) 


2.549 (0.169) 


— J.ZUZ 


u.uzo 


— 4. 141 


n m a 




Right 


2.361 (0.141) 


2.523 (0.224) 


2.610 (0.145) 


-3.914 


0.019 


-4.635 


0.009 


Inferior parietal lobes 


Left 


2.678 (0.175) 


2.826(0.186) 


2.922 (0.201) 


-2.845 


0.033 


-4.843 


0.007 




Right 


2.692 (0.177) 


2.728 (0.182) 


2.812 (0.159) 


-2.323 


0.041 


-2.684 


0.035 


Superior parietal lobes 


Left 


2.254 (0.144) 


2.318 (0.145) 


2.368 (0.137) 


-1.093 


NS 


-2.963 


0.033 




Right 


2.322 (0.109) 


2.421 (0.157) 


2.468 (0.127) 


— 1 .527 


NS 


— 3.307 


0.032 


Lateral occipital gyrus 


Left 


2.815 (0.284) 


2.874 (0.267) 


2.941 (0.297) 


-1.368 


NS 


-1.788 


NS 




Right 


2.554(0.285) 


2.649 (0.253) 


2.757 (0.313) 


-3.048 


0.032 


-2.534 


0.040 


Anterior cingulate cortex (ACC) 


Left 


2.721 (0.224) 


2.961 (0.219) 


3.080 (0.344) 


-3.755 


0.028 


-3.409 


0.029 




Right 


2.628 (0.335) 


2.882 (0.224) 


2.972 (0.348) 


-3.643 


0.021 


-4.083 


0.017 


Posterior cingulate cortex (PCC) 


Left 


2.601 (0.113) 


2.756 (0.141) 


2.771 (0.189) 


-3.286 


0.029 


-3.309 


0.031 




Right 


2.607 (0.157) 


2.715 (0.115) 


2.799 (0.168) 


-3.048 


0.032 


-3.443 


0.025 


Cortical thickness thickening (PTSD > 


control) 
















Calcarine gyrus 


Left 


2.536(0.301) 


2.517 (0.388) 


2.461 (0.253) 


2.058 


0.047 


4.426 


0.011 




Right 


2.586 (0.240) 


2.563 (0.333) 


2.517 (0.254) 


1.579 


NS 


1.692 


NS 



PTSD, post-traumatic stress disorder; SD, standard deviation. Results significant at false discovery rate threshold 0.05; NS: not significant. 



cerebral cortex, and cortical thinning in those regions may be associ- The decreased gray matter volume and density in different regions 

ated with the development of recent onset PTSD and PTSD symptom were found by the VBM method (Bonne et aL, 2008; Jatzko et aL, 
severity. 2006; Nardo et al., 2010; Zhang et al., 2011), but the VBM method 




Fig. 1. Cortical thickness difference between the PTSD patients and the normal controls. Region abbreviations from FreeSurfer model: SFG, superior frontal gyrus; HiG, hippocampal gyrus; 
1PL, inferior parietal lobes; ACC, anterior cingulate cortex; PCC, posterior cingulate cortex; SFG, superior frontal gyrus; RMFG, rostral middle frontal gyrus; LOG, lateral occipital gyrus. The T 
value color was overlaid on the reconstruction image, and the color bar was shown at the bottom. T values in voxelwise statistical maps were calculated by a general linear model of cortical 
thickness at each voxel covarying for age and gender. Indicated regions in yellow illustrate where the cortices of PTSD patients are thicker than that of the normal controls at corrected 
p < 0.05 (FDR corrected), while the indicated region in blue illustrates where the cortices of PTSD patients are thinner than that of the normal controls at corrected p < 0.05 (FDR 
corrected). (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 
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Left PTSD vs non-PTSD Right 




Fig. 2. The cortical thickness difference between PTSD group and non-PTSD group. Region abbreviations from FreeSurfer model: SFG, superior frontal gyrus; HiC, hippocampal gyrus; 1PL, 
inferior parietal lobes; ACC, anterior cingulate cortex; PCC, posterior cingulate cortex; SFG, superior frontal gyrus; RMFG, rostral middle frontal gyrus; LOG, lateral occipital gyrus. The T 
value color was overlaid on the reconstruction image, and the color bar was shown at the bottom. The age and gender effects were removed by regression. 



does not have the specificity for the highly convoluted structures of the 
cortex (Makris et al„ 2006). Surface-based morphometry is a highly 
reliable MRI-derived neuroanatomical measurement developed recently, 

Table 3 

Correlations of CAPS scores with regional cortical thickness. 



Regions 


Peak 

coordinates 
(X, Y, Z) 


Pearson 
r value 


P 

value 


Thickness at 
peak voxel 


Mean thickness of 
ROl (SD) mm 


Right hemisphere 










Inferior 


(-13,-79,18) 


-0.524 


0.007 


2.675 


2.717 (0.18) 


parietal 












lobes 












Superior 


(-10,91,-23) 


-0.484 


0.014 


2.384 


2.395 (0.21) 


frontal 












lobes 












Hippocampal 


(2,-24,-48) 


-0.422 


0.035 


2.289 


2.298 (0.13) 


gyrus 












ACC 


(35,73,-7) 


-0.473 


0.017 


2.725 


2.811 (0.31) 


PCC 


(31,-74,-14) 


-0.455 


0.024 


2.667 


2.713 (0.24) 


Left hemisphere 










Inferior 


(21,-65,18) 


-0.407 


0.043 


2.593 


2.613 (0.13) 


parietal 












lobes 












Superior 


(-6,51,47) 


-0.490 


0.013 


2.457 


2.499 (0.28) 


frontal 












lobes 












Hippocampal 


(0,-16,-52) 


-0.482 


0.015 


2.381 


2.393 (0.22) 


gyrus 












ACC 


(-34,41,21) 


-0.472 


0.017 


2.629 


2.631 (0.24) 


PCC 


(4,-64,-14) 


-0.413 


0.039 


2.651 


2.681 (0.19) 



CAPS, Clinician-Administered PTSD Scale. Results significant at false discovery rate 
threshold p < 0.05; PCC, posterior cingulate cortex; ACC, Anterior cingulate cortex; 
x, y, and z coordinates were located the peak vertex by Qdec surface coordinates. 



which overcomes the shortcomings of volume-based analysis, and is 
unaffected by voxel geometry (Wonderlick et al., 2009). The multiple- 
plane surface-based analysis keeps the advantage of surface-based analy- 
sis over volume-based analysis (De Leebeeck et al., 2007), and can accu- 
rately show the signs of cortical morphological compensatory plasticity 
in recent onset PTSD (Wang et al., 2010). 

The cortical thinning observed in the present study is partially consis- 
tent with the observations in previous studies about cortical thickness in 
recent onset or chronic PTSD (Geuze et al., 2008; Liu et al„ 2012; Rauch 
et al„ 2003; Woodward et al., 2009b), but is different from the observa- 
tions in some previous studies (Landre et al., 2010; Lyoo et al„ 2011). 
This may be due to the following reasons: first, in previous studies, the 
victim suffered from different traumas, such as combat related trauma 
(Geuze et al., 2008; Woodward et al., 2009b), sexual abuse (Landre 
et al., 201 0) and subway disaster (Lyoo et al„ 201 1 ). The different trauma 
types and trauma duration time may have different effects on the 
subjects' brain structures, and thus may cause the heterogeneity of 
study results. Second, in previous studies, the comparisons are be- 
tween PTSD and non PTSD (Geuze et al., 2008; Liu et al., 2012; 
Woodward et al„ 2009b), or between PTSD and the normal controls 
(Landre et al., 2010; Lyoo et al., 2011), which may neglect or mix 
the stress effects with PTSD effect on brain structure. There is 
evidence to show that gray matter volume of the limbic region has de- 
creased in traumatized non PTSD subjects (Ganzel et al., 2008) which pro- 
vide evidence that the stress is associated with structural changes of 
cerebral cortex. Third, most previous studies were about chronic PTSD 
(Geuze et al., 2008; Landre et al„ 2010; Lyoo et al., 2011; Woodward 
et al., 2009b). Their results reflected the pathogenesis of chronic PTSD, 
and the increased cortical thickness was found in recovery from chronic 
PTSD (Lyoo et al., 2011), possibly causing the inconsistency in the study 
results. In the present study, the subjects share the same traumatic 
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event, avoiding the above-mentioned confounding factors, and the cor- 
tical thickness was measured 6 months after trauma, so it is very likely 
to ascertain that the cortical thickness reductions in the current study 
were associated with recent onset PTSD caused by a single prolonged 
trauma. 

The hippocampus has been shown to play a role in contextual fear 
conditioning, and context modulation of extinction retention (Corcoran 
and Maren, 2001 ), and it is crucially involved in the encoding, consolida- 
tion and retrieval of declarative memories (Bekes, 2010). Dysfunction of 



the hippocampus is thought to be a central mechanism in the develop- 
ment of PTSD (Karl et al., 2006; Liberzon and Sripada, 2008; Wignall 
et al., 2004). Smaller bilateral hippocampal volume relative to the normal 
controls was also found in non PTSD subjects (Etldn and Wager, 2007), 
suggesting that the stress also has an effect on the hippocampus. Similar 
results were found in both recent onset PTSD (Zhang et al., 2011) and 
chronic PTSD (Bonne et al., 2008), and they were also verified in animal 
study (Golub et al., 201 1 ). In the current study, we found bilateral hippo- 
campus cortical thickness reductions in recent onset PTSD that is different 
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from previous studies (Geuze et al., 2008; Liu et al., 2012; Woodward 
et al., 2009b), and may be associated with severe stress or recent onset 
PTSD. The cortical thinning rate of the bilateral hippocampus in recent 
onset PTSD patients is much larger than that in non-PTSD subjects. 
These results imply that severe stress has profoundly selective effects on 
the hippocampus cortex, and may be associated with PTSD symptom 
severity. 

The present study shows that the bilateral cingulate cortex including 
ACC and PCC is thinning in recent PTSD, that is partially consistent with 
the results of meta-analyses of voxel-based morphometry studies of 
PTSD (Karl et al., 2006; Rauch et al., 2003). However, for the first time, 
we report cortical thinning in bilateral PCC in recent onset PTSD. The 
ACC plays a key role in frontolimbic networks associated with emotional 
regulating and cognitive functions (Paus, 2001). The cortical thickness 
of ACC correlates with PTSD symptom improvement (Dickie et al., 
2012), implying that the functional recovery of PTSD correlates 
with ACC cortical thickness. PCC is also related to a series of higher 
cognitive functions involved in PTSD, such as episodic and autobio- 
graphical memory, and self-referential processing (Liberzon and 
Sripada, 2008). Moreover, it has been suggested that PCC could be 
implicated in coping with physical threats, and plays a role in the 
processing of distressing material. Recently, fMRI studies have shown 
that spontaneous activity in the default network during rest is altered 
in PTSD (Bluhm et al., 2009; Wu et al., 2011). Thus, the PCC cortical 
thinning in the present study may suggest dysfunction in the processing 
of distressing material in severe trauma-exposed victims with and 
without recent onset PTSD. The average cortical thickness of the cingu- 
late cortex decreased almost 3 times in recent onset PTSD than in non- 
PTSD. This result implies that severe stress has profoundly selective 
effects on the bilateral cingulate cortex. 

In the present study, we find bilateral cortical thinning in the 
superior and inferior parietal and superior frontal lobe. The frontal 
lobe associates with the function of learning and memory, and the pari- 
etal lobe has extensive connections with all these brain regions 
(Liberzon and Sripada, 2008). The structural alteration observed in 
the inferior frontal lobe is related to fear inhibition, while the medial 
frontal cortex is related to conditional fear elimination (Jelicic and 
Merckelbach, 2004). Previous fMRI study also provided evidence 
that traumatic emotional reactions have a primary impact on the 
prefrontal systems engaged by the conscious processing of fear 
signals (Williams et al., 2006). Moreover, PTSD patients had abnor- 
mal frontal and parietal activity during working memory (Etkin 
and Wager, 2007), and a lower cerebral perfusion in the right 
parietal cortex (Chung et al., 2006; Schuff et al., 2011). Some authors 
suggest that apparent memory deficits in PTSD may due to impair- 
ment in attention and concentration in the frontal lobe, rather than 
in hippocampus dysfunction (Rossi et al., 2006). 

There was cortical thinning in the right lateral occipital gyrus, and 
thickening in the left calcarine cortex in recent onset PTSD. These 
regions are involved in visual information processing. It has been 
proved that the visual cortex has morphological compensatory plas- 
ticity in acquired strong visual stimulation (Noppeney et al., 2005). 
The fire disaster (Bronner et al., 2008; Chen et al., 2006) and mine- 
related trauma have strong visual stimulations (Zhang et al., 2011), 
which may also result in visual-cortex thinning. This result may associ- 
ate with the type of trauma. The survivors in the current study were 
trapped in a total darkness environment for 75 h, thus cortical thickness 
increase may reflect the compensatory effect of the visual cortex. In 
addition, lateral occipital gyrus may be involved in the higher-order cog- 
nitive functions (Cavanna and Trimble, 2006), the cortical thinning of 
lateral occipital gyrus may result cognitive impairment. 

We also find that CAPS is significantly correlated with the cortical 
thickness of the region and that the cortical thickness had significant 
changes, such as the inferior parietal and superior frontal lobes and 
the limbic system. Furthermore, the different regions have different 
scales of relevance. This is consistent with previous studies (Bryant 



et al., 2008; Liu et al., 2012), and may be related to the loss of neurons, 
as the duration of stress may be positively correlated with the degree 
of dendritic atrophy (Bonne et al., 2008). 

We acknowledge that the present study has several limitations that 
will be addressed in future study. First the sample size of the current 
study is small, thus statistical issues may obscure the subtle changes of 
cortical thickness. Second, non PTSD subjects in our study are all male, 
although the related study did not find a gender effect on hippocampus 
volume, especially between PTSD and trauma-exposed controls (Nardo 
et al., 2010). It would be better to include female subjects in the non 
PTSD group. Third, the average years in education of our study subjects 
are much lower than that of the previous studies. There is a report that 
lower pre-trauma intelligence increases the risk of developing PTSD 
symptoms (McFarlane, 1997) which may have an impact on the current 
study results. Moreover, we get the cortical thickness values from the 
older version of the software package (v4.0.2), which has shown that 
the calculations, given the same input, can vary up to 15% between 
different versions of the software in recent research (Gronenschild 
et al., 2012). In order to investigate the regionally subtle cortical thick- 
ness alterations, we would process the data by FreeSurfer version later 
than V5.2.0 in future research. 

In conclusion, the current study provides evidence that recent onset 
PTSD has a distinct pattern of cortical thinning in the limbic system, the 
frontal and parietal regions, and that may associate with pathophysiol- 
ogy of recent onset PTSD caused by a single prolonged trauma exposure. 
The severe stress still has an effect on the cerebral cortex in the same 
region as recent onset PTSD but at a different scale in non PTSD. Further- 
more, cortical alterations of limbic structural, parietal and frontal cortex 
correlate with current PTSD symptom severity. These findings suggest 
that the severe stress and recent onset PTSD pathogenesis is not a 
regional specific but a different scale in cortical thinning. The cortical 
thinning region is much larger than that found in VBM study, suggesting 
that the cortical thickness may be a more sensitive biomarker than gray 
matter volume in exploring pathophysiological change of recent onset 
PTSD. 
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